home *** CD-ROM | disk | FTP | other *** search
/ bioinformatics.org / bioinformatics.org_software.tar / www.bioinformatics.org / download / ecell2 / ecell220setup.exe / {app} / standard / STDR / RandomBiBiReactor.rd < prev    next >
Text File  |  2000-03-03  |  12KB  |  60 lines
  1. @CLASSNAME:RandomBiBiReactor
  2. @BASECLASS:FluxReactor
  3. @BRIEF_DESCRIPTION: Random Bi Bi reaction system
  4.  
  5. %VERSION: ecs-v09, 0.1
  6. @DATE: 30/6/1999
  7.  
  8. @AUTHOR: Masayuki Okayama
  9. @EMAIL: t98204mo@sfc.keio.ac.jp
  10.  
  11. @DESCRIPTION: A reactor class for a random bireactant system in the Bi Bi direction.
  12.  
  13. @EQUATION:
  14. {\tiny
  15. $v=\frac{M(P)\:[E]}{Den} :
  16. M(P)=k_{+7}k_{+8}(k_{+1}k_{-2}k_{+3}k_{+6}+k_{-1}k_{+2}k_{+4}k_{+6}+k_{+1}k_{-2}k_{+3}k_{+5}+k_{-1}k_{+2}k_{+4}k_{+5})[S_{1}][S2]+k_{-1}k_{-2}(-k_{-3}k_{-6}k_{+7}k_{-8}-k_{-4}k_{-6}k_{+7}k_{-8}-k_{-3}k_{-5}k_{-7}k_{+8}-k_{-4}k_{-5}k_{-7}k_{+8})[P1][Q]+k_{-5}k_{+6}k_{+8}(k_{+1}k_{-2}k_{+3}+k_{-1}k_{+2}k_{+4})[S_{1}][S2][Q]+k_{+1}k_{+3}k_{+4}k_{+7}k_{+8}(k_{+6}+k_{+5})[S_{1}]^{2}[S2]+k_{+2}k_{+3}k_{+4}k_{+7}k_{+8}(k_{+6}+k_{+5})[S_{1}][S2]^{2}+k_{-1}k_{-2}k_{-5}k_{-6}k_{-8}(-k_{-3}-k_{-4})[P1][Q]^{2}+k_{-1}k_{-3}k_{+4}(-k_{-6}k_{+7}k_{-8}-k_{-5}k_{-7}k_{+8})[S_{1}][P1][Q]+k_{-2}k_{+3}k_{-4}(-k_{-6}k_{+7}k_{-8}-k_{-5}k_{-7}k_{+8})[S2][P1][Q]+k_{+5}k_{-6}k_{+7}(k_{+1}k_{-2}k_{+3}+k_{-1}k_{+2}k_{+4})[S_{1}][S2][P1]+k_{-1}k_{-2}k_{-5}k_{-6}k_{-7}(-k_{-3}-k_{-4})[P1]^{2}[Q]+k_{+1}k_{+3}k_{+4}k_{-5}k_{+6}k_{+8}[S_{1}]^{2}[S2][Q]+k_{+2}k_{+3}k_{+4}k_{-5}k_{+6}k_{+8}[S_{1}][S2]^{2}[Q]-k_{-1}k_{-3}k_{+4}k_{-5}k_{-6}k_{-8}[S_{1}][P1][Q]^{2}-k_{-2}k_{+3}k_{-4}k_{-5}k_{-6}k_{-8}[S2][P1][Q]^{2}+k_{+1}k_{+3}k_{+4}k_{+5}k_{-6}k_{+7}[S_{1}]^{2}[S2][P1]+k_{+2}k_{+3}k_{+4}k_{+5}k_{-6}k_{+7}[S_{1}][S2]^{2}[P1]-k_{-1}k_{-3}k_{+4}k_{-5}k_{-6}k_{-7}[S_{1}][P1]^{2}[Q]-k_{-2}k_{+3}k_{-4}k_{-5}k_{-6}k_{-7}[S2][P1]^{2}[Q]$}
  17. {\tiny
  18. $Den=k_{-1}k_{-2}k_{+7}k_{+8}(k_{-3}+k_{-4}+k_{+5}+k_{+6})+k_{+7}k_{+8}(k_{+1}k_{-2}k_{-3}+k_{+1}k_{-2}k_{-4}+k_{+1}k_{-2}k_{+5}+k_{+1}k_{-2}k_{+6}+k_{-1}k_{-3}k_{+4}+k_{-1}k_{+4}k_{+5}+k_{-1}k_{+4}k_{+6})[S_{1}]+k_{-1}k_{-2}(k_{-3}k_{-5}k_{+8}+k_{-4}k_{-5}k_{+8}+k_{-5}k_{+6}k_{+8}+k_{-3}k_{+7}
  19. k_{-8}+k_{-4}k_{+7}k_{-8}+k_{+5}k_{+7}k_{-8}+k_{+6}k_{+7}k_{-8})[Q]+k_{-1}k_{-2}(k_{-3}k_{-6}k_{+7}+k_{-4}k_{-6}k_{+7}+k_{+5}k_{-6}k_{+7}+k_{-3}k_{-7}k_{+8}+k_{-4}k_{-7}k_{+8}+k_{+5}k_{-7}k_{+8}+k_{+6}k_{-7}k_{+8})[P1]+k_{+7}k_{+8}(k_{-1}k_{+2}k_{-3}+k_{-1}k_{+2}k_{-4}+k_{-1}k_{+2}k_{+5}+k_{-1}k_{+2}k_{+6}+k_{-2}k_{+3}k_{-4}+k_{-2}k_{+3}k_{+5}+k_{-2}k_{+3}k_{+6})[S2]+(k_{+1}k_{-2}k_{-3}k_{-5}k_{+8}+k_{+1}k_{-2}k_{-4}k_{-5}k_{+8}+k_{+1}k_{-2}k_{-5}k_{+6}k_{+8}+k_{-1}k_{-3}k_{+4}k_{-5}k_{+8}+k_{-1}k_{+4}k_{-5}k_{+6}k_{+8}+k_{-1}k_{-3}k_{+4}k_{+7}k_{-8}+k_{-1}k_{+4}k_{+5}k_{+7}k_{-8}+k_{-1}k_{+4}k_{+6}k_{+7}k_{-8})[S_{1}][Q]+(k_{+1}k_{-2}k_{-3}k_{-6}k_{+7}+k_{+1}k_{-2}k_{-4}k_{-6}k_{+7}+k_{+1}k_{-2}k_{+5}k_{-6}k_{+7}+k_{-1}k_{-3}k_{+4}k_{-6}k_{+7}+k_{-1}k_{+4}k_{+5}k_{-6}k_{+7}+k_{-1}k_{-3}k_{+4}k_{-7}k_{+8}+k_{-1}k_{+4}k_{+5}k_{-7}k_{+8}+k_{-1}k_{+4}k_{+6}k_{-7}k_{+8})[S_{1}][P1]+(k_{+1}k_{-2}k_{+3}k_{+5}k_{+8}+k_{+1}k_{-2}k_{+3}k_{+6}k_{+7}+k_{+1}k_{-2}k_{+3}k_{+7}k_{+8}+k_{+1}k_{+3}k_{-4}k_{+7}k_{+8}+k_{-1}k_{+2}k_{+4}k_{+5}k_{+8}+k_{-1}k_{+2}k_{+4}k_{+6}k_{+7}+k_{-1}k_{+2}k_{+4}k_{+7}k_{+8}+k_{+2}k_{-3}k_{+4}k_{+7}k_{+8}+k_{+3}k_{+4}k_{+5}k_{+7}k_{+8}+k_{+3}k_{+4}k_{+6}k_{+7}k_{+8})[S_{1}][S2]+k_{+1}k_{+4}k_{+7}k_{+8}(k_{-3}+k_{+5}+k_{+6})[S_{1}]^{2}+(k_{-1}k_{-2}k_{-3}k_{-5}k_{-6}+k_{-1}k_{-2}k_{-4}k_{-5}k_{-6}+k_{-1}k_{-2}k_{-5}k_{-7}k_{+8}+k_{-1}k_{-2}k_{-5}k_{+6}k_{-7}+k_{-1}k_{-4}k_{-5}k_{-7}k_{+8}+k_{-2}k_{-3}k_{-5}k_{-7}k_{+8}+k_{-1}k_{-2}k_{+5}k_{-6}k_{-8}+k_{-1}k_{-2}k_{-6}k_{+7}k_{-8}+k_{-1}k_{-4}k_{-6}k_{+7}k_{-8}+k_{-2}k_{-3}k_{-6}k_{+7}k_{-8})[P1][Q]+(k_{-1}k_{+2}k_{-3}k_{-5}k_{+8}+k_{-1}k_{+2}k_{-4}k_{-5}k_{+8}+k_{-1}k_{+2}k_{-5}k_{+6}k_{+8}+k_{-2}k_{+3}k_{-4}k_{-5}
  20. k_{+8}+k_{-2}k_{+3}k_{-5}k_{+6}k_{+8}+k_{-2}k_{+3}k_{-4}k_{+7}k_{-8}+k_{-2}k_{+3}k_{+5}k_{+7}k_{-8}+k_{-2}k_{+3}k_{+6}k_{+7}k_{-8})[S2][Q]+(k_{-1}k_{+2}k_{-3}k_{-6}k_{+7}+k_{-1}k_{+2}k_{-4}k_{-6}k_{+7}+k_{-1}k_{+2}k_{+5}k_{-6}k_{+7}+k_{-2}k_{+3}k_{-4}k_{-6}k_{+7}+k_{-2}k_{+3}k_{+5}k_{-6}k_{+7}+k_{-2}k_{+3}k_{-4}k_{-7}k_{+8}+k_{-2}k_{+3}k_{+5}k_{-7}k_{+8}+k_{-2}k_{+3}k_{+6}k_{-7}k_{+8})[S2][P1]+k_{+2}k_{+3}k_{+7}k_{+8}(k_{-4}+k_{+5}+k_{+6})[S2]^{2}+k_{-1}k_{-2}k_{-6}k_{-7}(k_{-3}+k_{-4}+k_{+5})[P1]^{2}+k_{-1}k_{-2}k_{-5}k_{-8}(k_{-3}+k_{-4}+k_{+6})[Q]^{2}+(k_{+1}k_{-2}k_{-3}k_{-5}k_{-6}+k_{+1}k_{-2}k_{-4}k_{-5}k_{-6}+k_{-1}k_{-3}k_{+4}k_{-5}k_{-6}+k_{-1}k_{+4}k_{-5}k_{-7}k_{+8}+k_{-1}k_{+4}k_{-5}k_{+6}k_{-7}+k_{-3}k_{+4}k_{-5}k_{-7}k_{+8}+k_{-1}k_{+4}k_{+5}k_{-6}k_{-8}+k_{-1}k_{+4}k_{-6}k_{+7}k_{-8}+k_{-3}k_{+4}k_{-6}k_{+7}k_{-8})[S_{1}][P1][Q]+(k_{+1}k_{-2}k_{+3}k_{+5}k_{-6}+k_{+1}k_{-2}k_{+3}k_{-6}k_{+7}+k_{+1}k_{+3}k_{-4}k_{-6}k_{+7}+k_{-1}k_{+2}k_{+4}k_{+5}k_{-6}+k_{-1}k_{+2}k_{+4}k_{-6}k_{+7}+k_{+2}k_{-3}k_{+4}k_{-6}k_{+7}+k_{+3}k_{+4}k_{+5}k_{-6}k_{+7}+k_{+3}k_{+4}k_{+5}k_{-7}k_{+8}+k_{+3}k_{+4}k_{+6}k_{-7}k_{+8})[S_{1}][S2][P1]+(k_{+1}k_{-2}k_{+3}k_{-5}k_{+8}+k_{+1}k_{-2}k_{+3}k_{-5}k_{+6}+k_{+1}k_{+3}k_{-4}k_{-5}k_{+8}+k_{-1}k_{+2}k_{+4}k_{-5}k_{+8}+k_{-1}k_{+2}k_{+4}k_{-5}k_{+6}+k_{+2}k_{-3}k_{+4}k_{-5}k_{+8}+k_{+3}k_{+4}k_{-5}k_{+6}k_{+8}+k_{+3}k_{+4}k_{+5}k_{+7}k_{-8}+k_{+3}k_{+4}k_{+6}k_{+7}k_{-8})[S_{1}][S2][Q]+k_{+1}k_{+3}k_{+4}(k_{+5}k_{+8}+k_{+6}k_{+7}+k_{+7}k_{+8})[S_{1}]^{2}[S2]+k_{+1}k_{+4}k_{-5}k_{+8}(k_{-3}+k_{+6})[S_{1}]^{2}[Q]+k_{+1}k_{+4}k_{-6}k_{+7}(k_{-3}+k_{+5})[S_{1}]^{2}[P1]+(k_{-1}k_{+2}k_{-3}k_{-5}k_{-6}+k_{-1}k_{+2}k_{-4}k_{-5}k_{-6}+k_{-2}k_{+3}
  21. k_{-4}k_{-5}k_{-6}+k_{-2}k_{+3}k_{-5}k_{-7}k_{+8}+k_{-2}k_{+3}k_{-5}k_{+6}k_{-7}+k_{+3}k_{-4}k_{-5}k_{-7}k_{+8}+k_{-2}k_{+3}k_{+5}k_{-6}k_{-8}+k_{-2}k_{+3}k_{-6}k_{+7}k_{-8}+k_{+3}k_{-4}k_{-6}k_{+7}k_{-8})[S2][P1][Q]+k_{+2}k_{+3}k_{+4}(k_{+5}k_{+8}+k_{+6}k_{+7}+k_{+7}k_{+8})[S_{1}][S2]^{2}+k_{+2}k_{+3}k_{-5}k_{+8}(k_{-4}+k_{+6})[S2]^{2}[Q]+k_{+2}k_{+3}k_{-6}k_{+7}(k_{-4}+k_{+5})[S2]^{2}[P1]+k_{-5}k_{-6}k_{-7}(k_{-1}k_{-2}+k_{-1}k_{-4}+k_{-2}k_{-3})[P1]^{2}[Q]+k_{-1}k_{+4}k_{-6}k_{-7}(k_{-3}+k_{+5})[S_{1}][P1]^{2}+k_{-2}k_{+3}k_{-6}k_{-7}(k_{-4}+k_{+5})[S2][P1]^{2}+k_{-5}k_{-6}k_{-8}(k_{-1}k_{-2}+k_{-1}k_{-4}+k_{-2}k_{-3})[P1][Q]^{2}+k_{-1}k_{+4}k_{-5}k_{-8}(k_{-3}+k_{+6})[S_{1}][Q]^{2}+k_{-2}k_{+3}k_{-5}k_{-8}(k_{-4}+k_{+6})[S2][Q]^{2}+(k_{+1}k_{-2}k_{+3}k_{-5}k_{-6}+k_{+1}k_{+3}k_{-4}k_{-5}k_{-6}+k_{-1}k_{+2}k_{+4}k_{-5}k_{-6}+k_{+2}k_{-3}k_{+4}k_{-5}k_{-6}+k_{+3}k_{+4}k_{-5}k_{-7}k_{+8}+k_{+3}k_{+4}k_{-5}k_{+6}k_{-7}+k_{+3}k_{+4}k_{+5}k_{-6}k_{-8}+k_{+3}k_{+4}k_{-6}k_{+7}k_{-8})[S_{1}][S2][P1][Q]+k_{+1}k_{+3}k_{+4}k_{-6}(k_{+5}+k_{+7})[S_{1}]^{2}[S2][P1]+k_{+1}k_{+3}k_{+4}k_{-5}(k_{+8}+k_{+6})[S_{1}]^{2}[S2][Q]+k_{+1}k_{-3}k_{+4}k_{-5}k_{-6}[S_{1}]^{2}[P1][Q]+k_{+2}k_{+3}k_{+4}k_{-6}(k_{+5}+k_{+7})[S_{1}][S2]^{2}[P1]+k_{+2}k_{+3}k_{+4}k_{-5}(k_{+8}+k_{+6})[S_{1}][S2]^{2}[Q]+k_{+2}k_{+3}k_{-4}k_{-5}k_{-6}[S2]^{2}[P1][Q]+k_{+4}k_{-5}k_{-6}k_{-7}(k_{-1}+k_{-3})[S_{1}][P1]^{2}[Q]+k_{+3}k_{-5}k_{-6}k_{-7}(k_{-2}+k_{-4})[S2][P1]^{2}[Q]+k_{+3}k_{+4}k_{+5}k_{-6}k_{-7}[S_{1}][S2][P1]^{2}+k_{+4}k_{-5}k_{-6}k_{-8}(k_{-1}+k_{-3})[S_{1}][P1][Q]^{2}+k_{+3}k_{-5}k_{-6}k_{-8}(k_{-2}+k_{-4})[S2][P1][Q]^{2}+k_{+3}k_{+4}k_{-5}k_{+6}k_{-8}[S_{1}][S2][Q]^{2}+k_{+1}k_{+3}k_{+4}k_{-5}k_{-6}[S_{1}]^{2}[S2][P1][Q]+k_{+2}k_{+3}k_{+4}k_{-5}k_{-6}[S_{1}][S2]^{2}[P1][Q]+k_{+3}k_{+4}k_{-5}k_{-6}k_{-7}[S_{1}][S2][P1]^{2}[Q]+k_{+3}k_{+4}k_{-5}k_{-6}k_{-8}[S_{1}][S2][P1][Q]^{2}
  22. $}
  23.  
  24. %SUBSTANCE: Substrate,2,2,
  25. %SUBSTANCE: Product,2,2,
  26. %SUBSTANCE: Catalyst,1,1,
  27.  
  28. %PARAMETER:k1, Float,
  29. %PARAMETER:k_1, Float,
  30. %PARAMETER:k2, Float,
  31. %PARAMETER:k_2, Float,
  32. %PARAMETER:k3, Float,
  33. %PARAMETER:k_3, Float,
  34. %PARAMETER:k4, Float,
  35. %PARAMETER:k_4, Float,
  36. %PARAMETER:k5, Float,
  37. %PARAMETER:k_5, Float,
  38. %PARAMETER:k6, Float,
  39. %PARAMETER:k_6, Float,
  40. %PARAMETER:k7, Float,
  41. %PARAMETER:k_7, Float,
  42. %PARAMETER:k8, Float,
  43. %PARAMETER:k_8, Float,
  44.  
  45. @REACT_FUNC:
  46.  
  47. Float velocity = catalyst(0)->quantity();
  48. Float S1= substrate(0)->concentration();
  49. Float S2= substrate(1)->concentration();
  50. Float P1= product(0)->concentration();
  51. Float P2= product(1)->concentration();
  52.  
  53. Float MP = k7*k8*(k1*k_2*k3*k6+k_1*k2*k4*k6+k1*k_2*k3*k5+k_1*k2*k4*k5)*S1*S2+k_1*k_2*(k_3*k_6*k7*k_8-k_4*k_6*k7*k_8-k_3*k_5*k_7*k8-k_4*k_5*k_7*k8)*P1*P2+k_5*k6*k8*(k1*k_2*k3+k_1*k2*k4)*S1*S2*P2+k1*k3*k4*k7*k8*(k6+k5)*pow(S1,2)*S2+k2*k3*k4*k7*k8*(k6+k5)*S1*pow(S2,2)+k_1*k_2*k_5*k_6*k_8*(k_3-k_4)*P1*pow(P2,2)+k_1*k_3*k4*(k_6*k7*k_8-k_5*k_7*k8)*S1*P1*P2+k_2*k3*k_4*(k_6*k7*k_8-k_5*k_7*k8)*S2*P1*P2+k5*k_6*k7*(k1*k_2*k3+k_1*k2*k4)*S1*S2*P1+k_1*k_2*k_5*k_6*k_7*(k_3-k_4)*pow(P1,2)*P2+k1*k3*k4*k_5*k6*k8*pow(S1,2)*S2*P2+k2*k3*k4*k_5*k6*k8*S1*pow(S2,2)*P2-k_1*k_3*k4*k_5*k_6*k_8*S1*P1*pow(P2,2)-k_2*k3*k_4*k_5*k_6*k_8*S2*P1*pow(P2,2)+k1*k3*k4*k5*k_6*k7*pow(S1,2)*S2*P1+k2*k3*k4*k5*k_6*k7*S1*pow(S2,2)*P1-k_1*k_3*k4*k_5*k_6*k_7*S1*pow(P1,2)*P2-k_2*k3*k_4*k_5*k_6*k_7*S2*pow(P1,2)*P2;
  54. Float Den= k_1*k_2*k7*k8*(k_3+k_4+k5+k6)+k7*k8*(k1*k_2*k_3+k1*k_2*k_4+k1*k_2*k5+k1*k_2*k6+k_1*k_3*k4+k_1*k4*k5+k_1*k4*k6)*S1+k_1*k_2*(k_3*k_5*k8+k_4*k_5*k8+k_5*k6*k8+k_3*k7*k_8+k_4*k7*k_8+k5*k7*k_8+k6*k7*k_8)*P2+k_1*k_2*(k_3*k_6*k7+k_4*k_6*k7+k5*k_6*k7+k_3*k_7*k8+k_4*k_7*k8+k5*k_7*k8+k6*k_7*k8)*P1+k7*k8*(k_1*k2*k_3+k_1*k2*k_4+k_1*k2*k5+k_1*k2*k6+k_2*k3*k_4+k_2*k3*k5+k_2*k3*k6)*S2+(k1*k_2*k_3*k_5*k8+k1*k_2*k_4*k_5*k8+k1*k_2*k_5*k6*k8+k_1*k_3*k4*k_5*k8+k_1*k4*k_5*k6*k8+k_1*k_3*k4*k7*k_8+k_1*k4*k5*k7*k_8+k_1*k4*k6*k7*k_8)*S1*P2+(k1*k_2*k_3*k_6*k7+k1*k_2*k_4*k_6*k7+k1*k_2*k5*k_6*k7+k_1*k_3*k4*k_6*k7+k_1*k4*k5*k_6*k7+k_1*k_3*k4*k_7*k8+k_1*k4*k5*k_7*k8+k_1*k4*k6*k_7*k8)*S1*P1+(k1*k_2*k3*k5*k8+k1*k_2*k3*k6*k7+k1*k_2*k3*k7*k8+k1*k3*k_4*k7*k8+k_1*k2*k4*k5*k8+k_1*k2*k4*k6*k7+k_1*k2*k4*k7*k8+k2*k_3*k4*k7*k8+k3*k4*k5*k7*k8+k3*k4*k6*k7*k8)*S1*S2+k1*k4*k7*k8*(k_3+k5+k6)*pow(S1,2)+(k_1*k_2*k_3*k_5*k_6+k_1*k_2*k_4*k_5*k_6+k_1*k_2*k_5*k_7*k8+k_1*k_2*k_5*k6*k_7+k_1*k_4*k_5*k_7*k8+k_2*k_3*k_5*k_7*k8+k_1*k_2*k5*k_6*k_8+k_1*k_2*k_6*k7*k_8+k_1*k_4*k_6*k7*k_8+k_2*k_3*k_6*k7*k_8)*P1*P2+(k_1*k2*k_3*k_5*k8+k_1*k2*k_4*k_5*k8+k_1*k2*k_5*k6*k8+k_2*k3*k_4*k_5*k8+k_2*k3*k_5*k6*k8+k_2*k3*k_4*k7*k_8+k_2*k3*k5*k7*k_8+k_2*k3*k6*k7*k_8)*S2*P2+(k_1*k2*k_3*k_6*k7+k_1*k2*k_4*k_6*k7+k_1*k2*k5*k_6*k7+k_2*k3*k_4*k_6*k7+k_2*k3*k5*k_6*k7+k_2*k3*k_4*k_7*k8+k_2*k3*k5*k_7*k8+k_2*k3*k6*k_7*k8)*S2*P1+k2*k3*k7*k8*(k_4+k5+k6)*pow(S2,2)+k_1*k_2*k_6*k_7*(k_3+k_4+k5)*pow(P1,2)+k_1*k_2*k_5*k_8*(k_3+k_4+k6)*pow(P2,2)+(k1*k_2*k_3*k_5*k_6+k1*k_2*k_4*k_5*k_6+k_1*k_3*k4*k_5*k_6+k_1*k4*k_5*k_7*k8+k_1*k4*k_5*k6*k_7+k_3*k4*k_5*k_7*k8+k_1*k4*k5*k_6*k_8+k_1*k4*k_6*k7*k_8+k_3*k4*k_6*k7*k_8)*S1*P1*P2+(k1*k_2*k3*k5*k_6+k1*k_2*k3*k_6*k7+k1*k3*k_4*k_6*k7+k_1*k2*k4*k5*k_6+k_1*k2*k4*k_6*k7+k2*k_3*k4*k_6*k7+k3*k4*k5*k_6*k7+k3*k4*k5*k_7*k8+k3*k4*k6*k_7*k8)*S1*S2*P1+(k1*k_2*k3*k_5*k8+k1*k_2*k3*k_5*k6+k1*k3*k_4*k_5*k8+k_1*k2*k4*k_5*k8+k_1*k2*k4*k_5*k6+k2*k_3*k4*k_5*k8+k3*k4*k_5*k6*k8+k3*k4*k5*k7*k_8+k3*k4*k6*k7*k_8)*S1*S2*P2+k1*k3*k4*(k5*k8+k6*k7+k7*k8)*pow(S1,2)*S2+k1*k4*k_5*k8*(k_3+k6)*pow(S1,2)*P2+k1*k4*k_6*k7*(k_3+k5)*pow(S1,2)*P1+(k_1*k2*k_3*k_5*k_6+k_1*k2*k_4*k_5*k_6+k_2*k3*k_4*k_5*k_6+k_2*k3*k_5*k_7*k8+k_2*k3*k_5*k6*k_7+k3*k_4*k_5*k_7*k8+k_2*k3*k5*k_6*k_8+k_2*k3*k_6*k7*k_8+k3*k_4*k_6*k7*k_8)*S2*P1*P2+k2*k3*k4*(k5*k8+k6*k7+k7*k8)*S1*pow(S2,2)+k2*k3*k_5*k8*(k_4+k6)*pow(S2,2)*P2+k2*k3*k_6*k7*(k_4+k5)*pow(S2,2)*P1+k_5*k_6*k_7*(k_1*k_2+k_1*k_4+k_2*k_3)*pow(P1,2)*P2+k_1*k4*k_6*k_7*(k_3+k5)*S1*pow(P1,2)+k_2*k3*k_6*k_7*(k_4+k5)*S2*pow(P1,2)+k_5*k_6*k_8*(k_1*k_2+k_1*k_4+k_2*k_3)*P1*pow(P2,2)+k_1*k4*k_5*k_8*(k_3+k6)*S1*pow(P2,2)+k_2*k3*k_5*k_8*(k_4+k6)*S2*pow(P2,2)+(k1*k_2*k3*k_5*k_6+k1*k3*k_4*k_5*k_6+k_1*k2*k4*k_5*k_6+k2*k_3*k4*k_5*k_6+k3*k4*k_5*k_7*k8+k3*k4*k_5*k6*k_7+k3*k4*k5*k_6*k_8+k3*k4*k_6*k7*k_8)*S1*S2*P1*P2+k1*k3*k4*k_6*(k5+k7)*pow(S1,2)*S2*P1+k1*k3*k4*k_5*(k8+k6)*pow(S1,2)*S2*P2+k1*k_3*k4*k_5*k_6*pow(S1,2)*P1*P2+k2*k3*k4*k_6*(k5+k7)*S1*pow(S2,2)*P1+k2*k3*k4*k_5*(k8+k6)*S1*pow(S2,2)*P2+k2*k3*k_4*k_5*k_6*pow(S2,2)*P1*P2+k4*k_5*k_6*k_7*(k_1+k_3)*S1*pow(P1,2)*P2+k3*k_5*k_6*k_7*(k_2+k_4)*S2*pow(P1,2)*P2+k3*k4*k5*k_6*k_7*S1*S2*pow(P1,2)+k4*k_5*k_6*k_8*(k_1+k_3)*S1*P1*pow(P2,2)+k3*k_5*k_6*k_8*(k_2+k_4)*S2*P1*pow(P2,2)+k3*k4*k_5*k6*k_8*S1*S2*pow(P2,2)+k1*k3*k4*k_5*k_6*pow(S1,2)*S2*P1*P2+k2*k3*k4*k_5*k_6*S1*pow(S2,2)*P1*P2+k3*k4*k_5*k_6*k_7*S1*S2*pow(P1,2)*P2+k3*k4*k_5*k_6*k_8*S1*S2*P1*pow(P2,2);
  55.  
  56. velocity *= MP;
  57. velocity /= Den;
  58.  
  59. process(velocity);
  60.